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In the first part of this paper, we study the spin-S Kitaev model using spin wave theory. We 
discover a remarkable geometry of the minimum energy surface in the TV-spin space. The classical 
ground states, called Cartesian or CN-ground states, whose number grows exponentially with the 
number of spins TV, form a set of points in the TV-spin space. These points are connected by a 
network of flat valleys in the TV-spin space, giving rise to a continuous family of classical ground 
states. Further, the CN-ground states have a correspondence with dimer coverings and with self 
avoiding walks on a honeycomb lattice. The zero point energy of our spin wave theory picks out a 
subset from a continuous family of classically degenerate states as the quantum ground states; the 
number of these states also grows exponentially with TV. In the second part, we present some exact 
results. For arbitrary spin-.?, we show that localized Z2 flux excitations are present by constructing 
plaquette operators with eigenvalues ±1 which commute with the Hamiltonian. This set of commut- 
ing plaquette operators leads to an exact vanishing of the spin-spin correlation functions, beyond 
nearest neighbor separation, found earlier for the spin- 1/2 model @. We introduce a generalized 
Jordan- Wigner transformation for the case of general spin-S, and find a complete set of commuting 
link operators, similar to the spin-1/2 model, thereby making the Z2 gauge structure more manifest. 
The Jordan- Wigner construction also leads, in a natural fashion, to Majorana fermion operators for 
half-odd-integer spin cases and hard-core boson operators for integer spin cases, strongly suggest- 
ing the presence of Majorana fermion and boson excitations in the respective low energy sectors. 
Finally, we present a modified Kitaev Hamiltonian which is exactly solvable for all half-odd-integer 
spins; it is equivalent to an exponentially large number of copies of spin-1/2 Kitaev Hamiltonians. 

PACS numbers: 75.10.Jm 



I. INTRODUCTION 

Frustrated quantum spin systems have become a new 
paradigm in condensed matter science. More and more 
systems are joining this family. A richness is emerging in 
terms of novel ground states and excitations. They not 
only enrich the basic science of strongly correlated elec- 
tron systems, but also have started playing a key role in 
an unexpected corner, namely quantum computers. Ki- 
taev has suggested that excitations of frustrated quantum 
spin systems have a special robustness, arising from their 
non-trivial topological property, which make them suit- 
able elements of a topological quantum computer. An ex- 
actly solvable two-dimensional frustrated spin-1/2 model 
introduced by Kitaev [l], Q exemplifies this. 

The spin-1/2 Kitaev model is interesting in its own 
right as a condensed matter spin model 0, 0,1E B IzL H, • 
In fact a similar model, called the compass model, al- 
though not exactly solvable, was introduced by Kugel 
and Khomskii in the late 1970s [k| to understand the 
magnetic properties of transition metal oxides which have 
orbital degeneracies. Recently an optical lattice realiza- 
tion of the spin-1/2 Kitaev model has been discussed [TT[ | . 

The spin- 5 Kitaev model for S > 1/2 is not ex- 
actly solvable. It is a challenging question if the Z% 
gauge structure and the presence of low energy Majorana 
fermions, discovered by Kitaev for the spin-1/2 case, sur- 
vive for arbitrary spin-S*. Are there differences between 
half-odd-integer and integer spins? In the present paper 



we approach the problem on two fronts. First we find ex- 
act classical ground states and perform a spin wave anal- 
ysis. The structure of the ground state manifold in the 
TV-spin space is rich: an exponentially large number of 
isolated points are connected by flat valleys. We will call 
the ground state corresponding to these isolated points as 
Cartesian or CN-ground states, as any given spin points 
along one of the three Cartesian directions. We also find 
the phenomenon of order from disorder in our spin wave 
analysis. We discover a nice connection between finding 
the CN-ground states and the dimer covering problem on 
the honeycomb lattice. In our spin wave analysis we find 
an equal number of finite frequency and zero frequency 
spin wave modes which live on self avoiding walks that 
are uniquely connected to the dimer covering. 

On the other front, we get some useful and exact results 
for the spin-S 1 Kitaev model, which prove the survival of 
the Z 2 gauge structure. There are also good indications 
that the low energy Majorana fermion excitations survive 
for half-odd-integer spins. Specifically we find plaquette 
operators (with eigenvalues ±1) which commute with the 
Hamiltonian. This set of commuting plaquette operators 
leads to a vanishing of the spin-spin correlation functions, 
beyond nearest neighbor separation, found earlier for the 
spin-1/2 model [fjj]. We also discover a Jordan- Wigner 
transformation for arbitrary S which leads to new bond 
operators (with eigenvalues ±1) which commute with the 
Hamiltonian. This makes the Z2 gauge structure mani- 
fest. The Jordan- Wigner construction leads, in a natural 
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fashion, to Majorana fermion operators for the case of 
half-odd-integer spins and to hard-core boson operators 
for the integer spin case, strongly suggesting the presence 
of Majorana fermions and bosons in the respective low 
energy sectors. 

The plan of the paper is as follows. In Sec. II, we 
consider a one-dimensional version of the Kitaev model 
and study spin wave theory in the large S limit. We 
find that although there is a continuous family of clas- 
sical ground states parametrized by an angle, the zero 
point energy of the spin waves picks out a discrete set of 
values of the angle as the quantum ground states; these 
correspond to dimer coverings on alternate bonds of the 
model. In Sec. Ill, we study the two-dimensional Ki- 
taev model using spin wave theory. We first present a 
general argument to find the classical ground state en- 
ergy. We then identify a discrete and infinite number 
of classical ground states; these have an interesting cor- 
respondence with self avoiding walks (SAW) and dimer 
coverings of the honeycomb lattice. The spin wave spec- 
trum is found to contain one- dimensional finite frequency 
spin wave modes and an equal number of zero frequency 
modes living on the SAWs. The zero point energy of the 
spin waves again picks out a special class of self avoid- 
ing walks as the quantum ground states. The number 
of classical and quantum ground states both grow expo- 
nentially with the number of sites, although the latter 
number grows slower than the former. In Sec. IV, wc 
construct, for any value of the spin S, an infinite set of Z 2 
operators which commute with each other and with the 
Hamiltonian. We use these operators to show that in any 
eigenstate of the Hamiltonian, the spin-spin correlations 
vanish unless the two spins are nearest neighbors; even 
for nearest neighbors, only certain components of the cor- 
relations are non-zero. Finally, we use a Jordan- Wigner- 
like transformation to construct a set of operators which 
act like Majorana fermions (hard-core bosons) for half- 
odd-integer (integer) values of S respectively. The Z 2 
operators defined earlier can be written as products of 
the Majorana fermion (hard-core boson) operators. In 
Sec. V, we present a modified Kitaev Hamiltonian whose 
energy spectrum can be found for any value of half-odd- 
integer spin; this model is equivalent to an exponentially 
large number of copies of the spin- 1/2 Kitaev Hamilto- 
nian. Some directions for future work are pointed out in 
Sec. VI. 



II. ONE-DIMENSIONAL KITAEV MODEL 

In this section, we will discuss a one-dimensional spin- 
s' model which is obtained by considering a single row of 
the Kitaev model in two dimensions. We illustrate the 
order from disorder phenomenon explicitly. The model 
is a spin S chain governed by the Hamiltonian 
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We assume that J > 0. (If J < 0, we can change its sign 
by performing a unitary rotation which flips the signs of 

Sf,i, s Zi> s !,2 and S h for a11 values of In E q- ©> the 

unit cells are labeled by i, and each unit cell has two spins 
labeled as 1 and 2. A factor of 1/S has been introduced 
in Eq. JT]) so that the ground state energy is proportional 
to S in the limit S — > 00. 

Let us introduce two vectors in the x — y plane, 



n = cos9 x + sin8 y 
and e = — sin9 x + cosO y. 



(2) 



Then a classical ground state of the Hamiltonian in Eq. 
(TT]) is given by the configuration 
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S h and Sf 2 



S n. 



(3) 



The classical energy of this state is E cl = -JSN/2, 
where N is the number of sites (the number of unit cells 
is N/2). Thus the classical ground states form a contin- 
uous family parametrized by an angle 9 which lies in the 
range [0, 2tt]. We will now perform a spin wave analysis 
and show that this picks out four values of 9 as having the 
lowest zero point energy; these correspond to Cartesian 
ground states. 

The spin wave spectrum around the ground state given 
in Eq. ((3]) can be found by using the Holstein-Primakoff 
(HP) transformation from spins to simple harmonic os- 
cillator raising and lowering operators [HI, [l3|, [ijj]. To 
obtain a HP Hamiltonian which is quadratic in bosons, 
we expand the fields as 

S*,i - Sn M - Pi ' 1 2 g g< ' 1 ^ + VS(egi,i + m,i) , 
Si, 2 = -Sh I 1 h2 nrt l ' 2 I + VS (-eqi,2 + zp h2 ) , 
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(4) 



where [qi, a ,Pj,b] — i$ijS a b- The spin wave Hamiltonian is 
then given by 



Hi, 



J J2 (Pli+Iti+Ph + I: 
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- cos 2 9 q it xq it 2 - sin 2 9 ^,2^+1,1) 
fc=0 v / \ 

where f(k) = - cos 2 9 - sin 2 9 e lk , (5) 
and k goes in steps of 4n/N. The spin wave energies are 



e k+ = Jv/TTL/m 

and e fc _ = Jy/l - \f(k)\, 



(6) 
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where |/(fc)| = yjl - sin 2 (28) sin 2 (k/2). The zero point 
energy is 

7T 

Ei, sw = jJ2(Vi+W)\ + vi - l/Wl) ■ (?) 

k=0 

We see that for each value of k, the spin wave energies 
have the same values for 9 and 7r/2 — 9. Now, 

- 5 M= - ~A=) < 0. (8) 

Thus the total spin wave energy ek++tk- increases mono- 
tonically as |/| decreases, i.e., as 9 increases from to 7r/4 
or decreases from tt/2 to 7r/4. Thus the zero point energy 
is minimum at 9 = 0, 7r/2, 7r and 3n/2, thereby picking 
out four points from the continuous family of classical 
ground states. These four points correspond to all the 
spins pointing along the ±x or ±y directions. 

Interestingly, the ground states chosen by the order 
from disorder phenomenon have two degenerate and non- 
dispersing spin wave branches with frequencies, €k+ = 
= J. It is easy to show that all these modes are 
localized on nearest neighbor bonds that have zero inter- 
action energy. 

Finally, let us briefly discuss the case in which the cou- 
plings are not equal on all the bonds. Suppose that the 
xx couplings have a strength J x and the yy couplings 
have a strength J y . If J x > J y , we find that the classical 
ground states are given by states in which the spins 1 and 
2 in each unit cell point in the +x, — x or — x, +x direc- 
tions. The classical ground state degeneracy is therefore 
2 N / 2 . We find that this de generacy is not broken by the 
zero point energy of the spin waves. 

III. TWO-DIMENSIONAL KITAEV MODEL 

We will now consider the spin-S* Kitaev model in two 
dimensions. This is a model on a honeycomb lattice with 
the Hamiltonian 

j-\-l—even 

+ S z 3l S z jl+1 ), (9) 

where j and I denote the column and row indices of the 
honeycomb lattice, respectively. We again assume, with- 
out loss of generality, that J > 0. Note that each spin 
is coupled to three other spins through xx, yy and zz 
couplings; we will denote the corresponding bonds as x, 
y and z bonds, respectively. (We present a schematic 
picture of the model in Fig. 1). We will first assume 
that the couplings on the three kinds of bonds are equal. 
Note that the honeycomb lattice is bipartite, with sites 



belonging to the two sublattices A and B having j + 1 as 
even and odd, respectively. If the total number of sites 
N is even, each sublattice has N/2 sites. 




FIG. 1: Schematic picture of the Kitaev model on a honey- 
comb lattice indicating the three kinds of bonds x, y and z. 
A hexagon with sites marked 1-6 is shown; the corresponding 
plaquette operator W p is defined in Eqs. (|16[) and (|18fl . 

We first present a general argument to obtain the clas- 
sical ground state energy of the Hamiltonian in Eq. 
We consider the spins at different sites S n to be classical 
(commuting) vectors, and introduce a Lagrange multi- 
plier A n at each site to enforce the relation S n = S 2 ; we 
do this by adding a term 

H * = - ^ E A » ( S n - ( 10 ) 
n 

to Eq. Q. Extremizing the sum of Eqs. ([9lfT0| leads to 
the equations 

oa \ na 

= -^n+m'S'n+m; (H) 

for any two neighboring sites n and n + m which are 
coupled by an a bond (a = x,y,z). Substituting Eq. 
(fTTj) in Eq. © and using the relation = S 2 , we find 
that the energy of such a state can be written in two ways 
which must be equal to each other, namely, 

E cl = JS A " = JS E V ( 12 ) 

n£A n£B 

Now, in any classical ground state, we can assume that 
for each site n, the spin on at least one of its three neigh- 
bors must point in such a direction that S^S^ +m 7^ 0. 
For such a pair, Eq. (TTTj) implies that A n A n + m = 1; note 
that n and n + m necessarily belong to different sublat- 
tices. Extending this argument to all pairs of neighboring 
sites, we conclude that A„ for all sites n belonging to sub- 
lattice A must have the same value, denoted as A^ , while 
A n for all sites n belonging to sublattice B must have 
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the same value, denoted as Xb, where XaXb = 1- Eq. 
(HU) then implies that E cl = (JSN/2)X A = (JSN/2)X B - 
The condition XaXb = 1 then implies that the minimum 
energy will be attained if Xa = Xb = — 1. Thus the clas- 
sical ground state is equal to —JSN/2 corresponding to 
A n = —1 at all sites. 

We will now explicitly find a large set of classical 
ground states. To this end, we observe the following 
interesting one to many correspondence of dimer cov- 
erings on a honeycomb lattice with a set of classical 
ground states of the Kitaev model that have identical en- 
ergy. Consider a covering of the honeycomb lattice with 
dimers, such that every site lies on a dimer. Associate a 
classical spin configuration to each dimer (bond) as fol- 
lows. Depending on whether it is an x, y or z bond, we 
put the two spins at the ends of the dimer as antiparal- 
lel and along the x, y or z direction respectively in spin 
space. This is the reason we call them Cartesian or CN- 
ground states. All these classical states have an identical 
energy — JSN/2. This follows from the fact that the in- 
teraction energy of the two spins of any dimer is —JS. 
Two neighboring spins not belonging to a dimer have zero 
interaction energy, either because they are orthogonal or 
the corresponding spin components do not appear in the 
bond interaction term. 

The number of dimer coverings on a honeycomb lattice 
has an asymptotic form (1,381) N / 2 fl5| . Further, the 
spins of each dimer can be in two possible antiparallcl 
states; hence we have 2 N I 2 classical spin configurations 
for each dimer covering. This makes the total degeneracy 
of the CN-ground states to be (1.662)^. 

An important question is whether these classical states 
remain stable under quantum fluctuations. We answer 
this question in two steps. First we show that the discrete 
and exponentially large set of degenerate states found 
above are further connected by flat valleys in the A-spin 
space. Then we perform a spin wave analysis and show 
that to the leading order in l/S there are no negative 
energy spin wave excitations. This ensures local stability 
of our quantum ground states. 

The discrete set of degenerate states obtained from the 
dimer coverings forms a set of isolated points in the N- 
spin space, i.e., {S 2 } N , where S 2 denotes the surface of a 
sphere in three dimensions. We will now show that there 
are flat valleys defined by a set of continuous parameter 
which connect the discrete points. To see this consider a 
set of self avoiding walks (SAWs) that completely covers 
the lattice, such that each lattice site appears on one and 
only one SAW. In each SAW, let alternate bonds form 
dimers. Each SAW must be either infinitely long or must 
be a closed loop consisting of an even number of bonds 
(this is because each bond on the honeycomb lattice goes 
from a site on sublattice A to a site on sublattice B). 

Let us now consider the sites and bonds lying on one 
particular SAW. To be specific, let us suppose that some- 
where in the middle of the SAW, we have some sites and 
bonds of the form ■ ■ ■ 1 — x — 2 — z — 3 — y — 4 — x — 5 ■ ■ ■ . 
A discrete classical ground state is then given by one 



in which the spins at the sites 1,2,3,4,5 point along 
z, — z, z, — x, x; this has an energy of —2JS for the four 
bonds xzyx. Another discrete classical ground state 
is given by taking the same five spins to point along 
x, — x, y,— y,y. We now note that a continuous family 
of classical ground states which interpolates between the 
above two discrete states is given by a configuration in 
which the five spins point along sin #x + cos 9i, — sin 0x — 
cos#z,sin(9y + cos#z, — cos#x — sinf^cos^x + sin(?y, 
where 9 goes from to 7r/2; the energy of the four bonds 
for this configuration is —2JS for all values of 9. This 
transformation can be extended to all the sites of the 
SAW. 

Thus we have a continuous transformation taking us 
from one discrete classical ground state of a SAW to an- 
other; we will call such a transformation a slide. A slide 
is parametrized by an angle 9, and it takes us from a dis- 
crete classical ground state in which the even numbered 
bond energies are minimized to one in which the odd 
numbered bond energies are minimized. We now observe 
that if a slide rotates the spin on a particular site in the 
x - y plane (such as site 4 in the previous paragraph), 
then that site must be coupled to a site in the neigh- 
boring SAW by a z bond; hence the classical interaction 
energy of that site to its neighboring site on the other 
SAW remains zero throughout the slide. This is true for 
any two neighboring sites belonging to different SAWs, 
no matter in which plane each of them is rotated during 
the slides of the two SAWs. Thus, a slide can be carried 
out on each SAW separately without changing the clas- 
sical interaction energy between the two SAWs; hence 
there is a continuous family of classical ground states on 
each SAW. However, based on the results in Sec. II, we 
expect that the zero point energy of the spin waves about 
such a continuous family of ground states in a SAW will 
be minimized for a discrete set of values of 9 which cor- 
responds to the spin at each site pointing along one of 
the six directions ±x, ±y or ±z. We will therefore con- 
sider below only the discrete set of classical ground states 
described in the previous sentence. 

We have seen that the discrete classical ground states 
correspond to dimer coverings of the honeycomb lattice. 
An interesting question to ask is whether all dimer cover- 
ings can be continuously connected to each other through 
the continuous families of classical ground states. We will 
prove that the answer is yes, by showing that any discrete 
classical ground state can be transformed by a succession 
of slides to a classical ground state in which all the dimers 
lie on the z bonds (we will call these vertical dimers) . For 
any dimer covering, the lattice can be covered by SAWs. 
We now choose these SAWs as follows. If all the dimers 
are vertical, there is nothing more to be done. If at least 
one dimer is non-vertical, we consider that dimer; each 
end of it also belongs to a vertical bond which is not a 
dimer (since no point can belong to two dimers). We go 
to the other end of that vertical bond; that end must 
belong to a non-vertical dimer. Continuing in this way, 
we get a SAW that consists of alternating non-vertical 
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dimers and vertical bonds which are not dimers. We now 
apply a slide to this SAW; we then get a SAW all of 
whose dimers are vertical. (Figure 2 shows an example 
of two SAWs which are connected by a slide). We then 
repeat the process of taking another non-vertical dimcr 
(which has no points in common with the previous SAW) , 
constructing a new SAW from it using the above proce- 
dure, and finally performing a slide which converts all the 
dimers to vertical ones. By repeating this until we have 
SAWs covering all the sites, we reach the state in which 
all the dimers are vertical. Thus all dimer coverings of 
the honeycomb lattice are connected through continuous 
families of classical ground states to the purely vertical 
dimer covering. Therefore all dimer coverings are also 
connected to each other through continuous families. 




FIG. 2: Pictures of two self avoiding walks, marked as 1, 
which are related to each other by a slide. The walk in (a) 
has only non-vertical dimers, while the walk in (b) has only 
vertical dimers. The dimers are shown by solid lines. 

We will now use the HP transformation to compute the 
spin wave spectrum about any one of the discrete classi- 
cal ground states given by a dimer covering. Consider a 
vertical dimer on which the two spins point along the ±z 
direction. Let us consider these two spins separately. 

(i) For the spin pointing along the z direction, 
we have S z = S - (p 2 + q 2 )/2, but (S x ,Sy) 
can be chosen in four different ways, namely, 
VS(q,p),VS(-q,-p),y/S(p,~q), and y/S(~p,q), up to 
the lowest order in the HP transformation. 

(ii) For the spin pointing along the — z direction, we have 
S z = -S + {p 2 + q 2 )/2, but (S x , S y ) can be chosen in 
four ways, namely, \/~S(p, q), \/S(~p, —q), VS(q, —p) and 
VS(-q,p), 

The term coupling the two spins at the opposite ends 
of the vertical dimer is given by SfSj = [S — (p 2 + 



q 2 )/2] [-S+(p 2 + q 2 )/2] ^-S 2 + (S/2)(p 2 + q 2 +p 2 + q 2 ) 
up to order S. Thus there is no coupling between the two 
spins to this order in S. 

I I I I I 

'1 

I I I I I 




FIG. 3: Picture of a set of self avoiding walks (dotted lines); 
one of the walks is marked as 1. The dimers are shown by 
solid lines. 

For any dimer covering of the honeycomb lattice, there 
is a set of self avoiding walks (SAW) covering the lattice 
such that none of the bonds appearing on a SAW is a 
dimer. Hence the two spins appearing at the two ends 
of any dimer belong to different SAWs, as shown in Fig. 
3. (Notice that these SAWs are different from the ones 
we discussed earlier while finding our flat valley through 
the slide operation; in those SAWs, alternate bonds were 
dimers). As we saw above, sites belonging to different 
SAWs are decoupled from each other up to order S. We 
can therefore carry out a spin wave analysis for each SAW 
separately. Now consider a SAW which forms a closed 
loop with n sites, where we saw earlier that n must be 
an even integer. (The minimum value of n is 6 corre- 
sponding to a hexagon). Let n — 2m. As we go around 
the loop, we choose the spin variables along the loops 
to be q and p alternately, so that the couplings between 
nearest neighbors involve either q m , q n or p m ,p n , but not 
q m ,Pn- Because of the two cases (i) and (ii) discussed 
above, the loop may have either periodic boundary con- 
dition (PBC) or antiperiodic boundary condition (ABC). 
Ignoring a constant, we find the spin wave Hamiltonian 
for the SAW to be 

T n 

h sw = ^ £ (p? + a 

i=l 
m — 1 

+ J X] (P2i-lP2i + fefe+l) 
»=1 

+ J (Pn-iPn ± q n Qx), (13) 

with either PBC or ABC for the last bond connecting 
sites n and 1; the sign of the q n qi term is + and — 
in the two cases respectively. The SAW has m unit 
cells, each consisting of 2 sites. The normal modes 
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can be characterized by a momentum k, where k = 
0, 27r/m, • • • , (27rm — 27r)/m in the case of PBC, and 
k = Tt/m, 3n/m, ••• , (27rm — ir)/m in the case of ABC. 
We now find the normal mode frequencies by solving the 
classical Hamiltonian equations of motion. For each mo- 
mentum k, we find that there are two frequencies given 
by and uJk = 2J\ cos(fc/2)|. The existence of a zero 
energy mode for each A: is a signature of the enormous 
ground state degeneracy at the classical level; the zero 
mode also implies that the spin wave correction to the 
expectation value of the spin at each site diverges. 

The spin wave normal modes are interesting. Even 
though they are formally characterized by a wave vec- 
tor, they are defined on self avoiding strings of varying 
shapes and sizes on the lattice. Secondly, by construc- 
tion, the classical interaction energy of the spins on a 
SAW is identically zero, bond by bond; it is neither a 
minimum nor a maximum. We have a collection of one- 
dimensional spin waves living on SAWs. Further these 
spin waves have a linear dispersion, uj-x+q ~ 2J|g| at low 
frequencies, around k — tt. 

A linear dispersion is known in the spin-1/2 Heisenberg 
antifcrromagnetic system in one dimension. However, the 
spin waves there are spin-1 excitations. In the Kitaev 
model under consideration, this linear dispersion occurs 
for both antiferromagnetic and ferromagnetic couplings 
and is non-degenerate, indicating some kind of spin-zero 
or real scalar character of the spin wave quanta. The frus- 
tration in the Kitaev model seems to induce an effective 
antiferromagnetic behavior along the SAW lines what- 
ever be the sign of J. This linear spin wave spectrum 
should be considered as a precursor to the linear Ma- 
jorana spectrum that one gets for the spin-1/2 Kitaev 
model. It is likely that these scalar spin wave quanta 
undergo quantum number fractionization leading to Ma- 
jorana fermions. 

We now calculate the zero point energy per site eo = 
J2k w fe/(2n) as a function of n > 6. We find that 

(i) for m odd, eo/J = cosec(ir / n) / n for PBC and 
cot(7r/n)/n for ABC. 

(ii) for m even, eo/J — cot(-7r/n)/n for PBC and 
cosec(ir / n) / n for ABC. 

In all cases, eo/J — * 1/ir ~ 0.318 as n — > oo. We find 
that the minimum value of eo occurs if n = 6 and we 
have ABC. In that case, e /J = V3/6 ~ 0.289. 

Now we show that the above minimum for the zero 
point energy per site can actually be achieved for the 
entire honeycomb lattice. Consider a dimer covering 
such that the corresponding SAWs cover the lattice with 
hexagons, as shown in Fig. 4. (There are three such 
dimer coverings; this will contribute a factor of 3 when 
we compute the ground state degeneracy below). For 
each such dimer covering, we take the two spins on all 
dimers on bonds of type n to point in the same way, say, 
the top spin pointing along n and the bottom spin point- 
ing along — n (here n = x,y or z). Then we find that 
each hexagon has ABC and therefore eo = 0.289J. 

Next, we note that for each of the three dimer cov- 




FIG. 4: Picture of a ground state in which the self avoiding 
walks (dotted lines) form hexagons; three of these are marked 
as 1-3. The dimers are shown by solid lines. 



erings, the SAWs cover the lattice with N/6 hexagons. 
These hexagons form a triangular lattice containing N/3 
triangles. We then see that the ABC of the Hamiltonian 
in Eq. (fT5|) continues to hold if we flip both the spins on 
all the three dimers which pairwise join three hexagons 
which form a triangle. So even after the dimer covering 
is fixed, we still have an exponentially large number of 
ground states given by 2 N / 3 corresponding to the differ- 
ent ways in which pairs of spins on different dimers can 
point. Since there are three possible dimer coverings, the 
total degeneracy of the quantum ground states is 3 x 

2 N/3 

which goes as (1.260)^ for large N. This is a smaller ex- 
ponential than the number of discrete classical ground 
states which goes as (1.662)^ as stated earlier. 

To conclude, we have shown that the dimer coverings 
which have the minimum zero point energy per site given 
by 0.289 J are the ones which correspond to SAWs cover- 
ing the lattice with hexagons. The number of quantum 
ground states picked out by the zero point energy of the 
spin waves still grows exponentially with the number of 
sites, but it grows more slowly than the number of clas- 
sical ground states. 

Finally, we can again consider what happens if the cou- 
plings on the three kinds of bonds are different, say, J x , 
J y and J z . We again find that if one of these, say J z , 
is larger than the other two, then the classical ground 
states are given by the state with purely vertical dimers; 
in each such dimer, the two spins can point along the 
+z, — z or — z, +z directions. Hence the number of clas- 
sical ground states is 1 N I 2 . We find that this degeneracy 
is not broken by the zero point energy of the spin waves. 



IV. CONSERVED Z 2 FLUXES 

In this section we construct commuting operators for 
the spin-S" Kitaev model and also generalize the Jordan- 
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Wigner transformation for the spin-,!? case. In the process 
we get exact results and some new insights. 

For the spin-1/2 Kitaev model in two dimensions, it is 
known that there is a conserved quantum number associ- 
ated with each hexagon. When the model is rewritten as 
a Zi gauge theory, these conserved quantities correspond 
to the Zi fluxes passing through the hexagon. Since the 
number of such quantum numbers (or hexagons) is N/2, 
the Hilbert space decomposes into 2 N I 2 independent sec- 
tors corresponding to each flux independently taking the 
values ±1. We will now show that all this continues to 
hold for arbitrary values of the spin S, integer or half- 
odd-integer. 

We first note that the three spin operators S x ,S y ,S z 
satisfy the identities 

e" s ° S b e~™ sa = S b if a = b, 

= - S b if a £ b, (14) 

for a,b = x,y,z. Now consider a hexagon h with sites 
labeled 1, • • • , 6, with the Hamiltonian 

TT _ QX QX | QV OV i QZ QZ , QX QX , Q& Qy I QZ QZ 

i±h — u^U2 t ^2 3 ' ^3^4 ' ^4^5 ' J t G ' 6 1" 

(15) 

If we define an operator 

W p = e i7r(s i + s * + s 3 + s| + st + si) ( 16 ) 

(see Fig. 1), then it follows from Eq. (fl4|) that 

WpH^Wp)- 1 = H h . (17) 

Thus W p commutes with the Hamiltonian H^. It is easy 
to check that W p also commutes with the other terms of 
the full Hamiltonian coming from other plaquettes. For 
S = 1/2, we observe that 

Wp = - a\ of of of of of, (18) 

where a a denote the Pauli matrices. 

Since e ±l27rSa = (-1) 2S , we have e Z7rSa = 
(-l) 2S e - i7rS "\ It then follows from Eq. (14]) that 

e i«S" e iTrS> = ^s^^S" if a ^ b (1Q) 

It then follows that W p will commute with W v < since the 
two will share an even number of sites. It also follows that 
Wp = 1. Thus the Wp's are a set of mutually commuting 
conserved operators with eigenvalues equal to ±1. 

We note that for S = 1/2, these operators are the same 
as the conserved flux operators. We will henceforth refer 
to these operators as the flux operators for all S. While 
we do not have a gauge theoretic formalism of the model 
for S > 1/2, we note that we can associate a conserved 
Z2 quantum number with every closed loop on the lattice 
as follows. For every site on the closed loop, we define 
the normal direction as the direction associated with the 
bond that does not belong to the loop. If (ii, i%, . . . , ijsr) 



are the sites for that loop and (ai, . . . , ajv) are the cor- 
responding normal directions, then the conserved quan- 
tity is 

N 

W L = H e i7rS % . (20) 

n=l 

It is interesting to note that these closed string operators 
are very similar to those defined by Den Nijs and Rom- 
melse [l6| in the context of S = 1 chains and generalized 
to arbitrary values of S by Oshikawa [TTJ • 

A. The ground state fluxes 

In the S = 1/2 model, it has been proven that the 
values of the flux operators are the same for all the ele- 
mentary hexagonal plaquettes and are equal to +1. Our 
semiclassical results indicate that this may be true for all 
values of S. 

Consider the quantum state constructed by taking the 
direct product of the spin coherent states corresponding 
to any of the classical ground state configurations de- 
fined by a dimer covering of the lattice. We will refer to 
these states as the semiclassical ground states. Consider 
a dimer covering which defines a set of closed self avoiding 
loops such that the dimers form the normal directions to 
these loops. The semiclassical ground state correspond- 
ing to such a dimer covering is a simultaneous eigenstate 
of the flux operators corresponding to the closed loops 
defined by the dimer covering. The eigenvalues are given 
by 

N 

II e " S? "" IV>> = e- 5 ^- (21) 

where p n = ±1 depending on the polarization of the spin 
at the n th site. Since the honeycomb lattice is bipar- 
tite, all closed loops have an even number of sites. The 
eigenvalues are thus e l27r ' s , where I is an integer. If S is 
an integer, then we see that the eigenvalues are always 
1. If S is a half-odd-integer, then it is (—1)'. It can be 
checked that the classical configurations which lead to 
ABC in the spin wave Hamiltonian correspond to even 

Thus we see that the spin wave fluctuations pick out 
the states which are simultaneous eigenstates of the flux 
operators corresponding to one third of the elementary 
plaquettes. It further picks out the states which have 
eigenvalue +1 for these operators. We note that the 
semiclassical ground states can never be the simultane- 
ous eigenstate of all the flux operators. The spin wave 
fluctuations seem to pick out the states which are the 
simultaneous eigenstates of the maximal number of ele- 
mentary plaquette flux operators with eigenvalue +1. In 
the S = 1/2 case, the exact ground state is a simulta- 
neous eigenstate of all the flux operators with eigenvalue 
+1. Our results indicate that this may be the case for all 
S also. 
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B. The flux basis and spin-spin correlations 

One of the intriguing features of the Kitaev model is 
the peculiar form of the spin-spin correlation functions 
@. Only the nearest neighbor correlations are non-zero. 
Further, only SfSJ on the x bonds, SfSj on the y bonds 
and Sf S? on the z bonds have non-zero values. These 
results are true for not only the ground state but also for 
any eigenstate of the model. In fact, as we will show, this 
form of the correlation function characterizes the set of 
simultaneous eigenstates of the flux operators. 

Consider any simultaneous eigenstate of all the elemen- 
tary flux operators denoted by where p n = ±1 is 
the eigenvalue of the flux operator of the n th plaquette. 
The spin operator Sf acting on this state produces an- 
other simultaneous eigenstate of the flux operators, with 
the eigenvalues of the two plaquettes which share the 
bond + a) flipped. This follows from Eqs. (fTO|) and 
(fT6|) . As two states with different sets of values of p n are 
orthogonal, the only non-zero spin-spin correlations are 

(S*S? +i ), (SfSf +i ), and <S?S^>. (22) 

Thus, as stated earlier, the flux operators define a basis 
of a peculiar kind of spin liquid. This basis is formally 
constructed in the fermionic formalism of the spin-1/2 
model. We will now show that it is possible to have an 
analogous Jordan- Wigner construction to construct this 
basis in the spin-5 1 model also. However, unlike the spin- 
1/2 case our construction will not exhaust all the states 
in the Hilbert space and thus does not lead to an exact 
solution of the model. 

Consider a Hamilton path running through the lattice. 
For simplicity, we consider an infinite lattice and take 
the path to go through the x and y bonds only. The 
construction can be generalized to any SAW [18| . Define 
the "disorder operators" on the n th site to be 

Min = Y[ e^ s L+S) , (23) 

where the 's commute with each other and = 1 . 
At each site i n , the Hamilton path is used to classify 
the three bonds as incoming, outgoing and normal. The 
normal bond for our path is the z bond. We denote 
the spin operator in the incoming and outgoing bond 
directions by S* 1 and 5' 2 respectively. We now define 
two operators at each site, 

6. = e M < +S W, 



It then follows that 

WinCim ( "0 £im£*« — $rnni 

Zi n Xi m - (-1) 2S X^&„ = 0. (25) 



Thus the £ and \ operators are Majorana fermions for 
half-odd-integer spins and hard-core bosons for integer 
spins. 

We now consider the commutators of the £j n and \i n 
with the Hamiltonian. Consider the terms £i„i?£i„ and 
Xi n HXi n - All t ne terms in the Hamiltonian involving 
spins at sites i m for m < n are left invariant. This is 
because either both or none of the spins have their signs 
flipped. So we have 

H£ in = H - 2 J t3 S\l S£_ x - 2 J a „ S£ S> +&n , 
Xi n H X i n = H-2J an S« : S« :+&n . (26) 

It then follows that the operators defined on the normal 
bonds, 

uy = ^ S XiXj, (27) 

form a set of mutually commuting operators which also 
commute with the Hamiltonian. The flux operator on 
any elementary plaquette is equal to the product of the 
Uij operators on the two z bonds of the hexagon just as 
in the 5=1/2 case. 

Thus, just as in the S = 1/2 case, the flux basis is 
easy to construct in terms of the Xi operators. However, 
writing it in a simple form in terms of the original spins 
remains a challenge. 

V. AN EXACTLY SOLVABLE HIGHER SPIN 
MODEL WITH FREE MAJORANA FERMIONS 

For the case of half-odd-integer spin S, we find that it 
is not possible to write the spin-S 1 Kitaev model in Eq. 
<j9j) in terms of local interactions between the Majorana 
fermions introduced above. The difficulty arises from not 
being able to invert Eq. (f!M]) to obtain the spin opera- 
tors. However, for any finite value of half-odd-integer 
spin, we can define a modified Kitaev Hamiltonian, for 
which we can give the exact spectrum and degeneracies. 
Our modified Kitaev Hamiltonian has the form 

h = | E &w + T UAi 

j-\-t=even 

+ ^j+i), (28) 

where r a = e l7lSa is the 7r-rotation operator introduced 
in Eq. IT4|) . Now, the operators r a only connect states 
which have the same magnitude of S z , in the S z basis; 
this is because e™ s * and e mSV acting on \S Z = m) give 
\S Z = —to). Hence the r a 's reduce to 2 x 2 blocks in the 
basis of eigenstates of S z ; the number of such blocks is 
equal to 5-1-1/2 corresponding to m — 1/2, 3/2, • ■ • , S. 
For example, in the case of S = 3/2, we find that in the 
basis of the two states with S z = ±3/2, t x = —io x , 
t v = ia y , and t z = —i<r z , while in the basis of the two 
states with S z = ±1/2, t x = -ia x , r» = -icr v , and 
t z = ia z . 
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In view of this, the (2S+1) N dimensional Hilbert space 
for N spins decomposes into (S + 1/2) N copies of 2 N di- 
mensional Hilbert spaces. Inside each copy, the Hamilto- 
nian in Eq. (f2"5| behaves exactly like the Kitaev Hamil- 
tonian in Eq. ([9]), leading to the identical spectrum and 
physical properties. The degeneracy of (S + 1/2) N of 
the Hamiltonian (which is not a gauge degeneracy) is an 
unusual decomposition of the eigenstates and may have 
some special use in quantum computation. 

In the case of integer spins, Eq. (flU)) shows that the 
r a operators commute with each other for all the values 
of a = x, y, z. Hence, they can be diagonalized simulta- 
neously; each of the diagonal entries is equal to ±1. The 
integer spin version of Eq. (|28|) therefore reduces to a 
kind of Ising spin model which is classical rather than 
quantum mechanical. 

VI. DISCUSSION 

In this paper we have presented a large spin analysis 
of the Kitaev model, whose spin-1/2 end is exactly solv- 
able. We find a classical ground state structure which has 
a non-trivial geometry in the TV-spin space. There are 
discrete sets of points that are connected by flat valleys. 



Our spin wave analysis gives either zero or positive energy 
excitations, indicating local stability of the degenerate 
set of vacuua. Further, depending on the vacuum cho- 
sen (which depends on the dimer covering pattern), the 
spin waves are localized on the SAW curves. It would be 
interesting to see how and when quantum number frac- 
tionization occurs and Majorana fermions emerge, when 
one goes beyond harmonic spin wave theory. 

A class of excitations that we have not studied in this 
paper is the one involving dimer coverings containing one 
or more defects. By definition, a defect site does not lie 
on a dimer. We can consider ground states about these 
defective dimer coverings; these will define ground states 
containing a topological spin defect. As the energy of one 
topological defect is of the order of JS, there is a finite 
gap to these excitations. It would be interesting to study 
whether such a defect might be related to a Majorana 
fermion excitation. 
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